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The Einstein relation for a driven moderately dense granular gas in d-dimensions is analyzed in 
the context of the Enskog kinetic equation. The Enskog equation neglects velocity correlations but 
retains spatial correlations arising from volume exclusion effects. As expected, there is a breakdown 
of the Einstein relation e = D/(To/-l) 7^ 1 relating diffusion D and mobility n, To being the tem- 
, perature of the impurity. The kinetic theory results also show that the violation of the Einstein 

relation is only due to the strong non-Maxwellian behavior of the reference state of the impurity 
particles. The deviation of e from unity becomes more significant as the solid volume fraction and 
the inelasticity increase, especially when the system is driven by the action of a Gaussian thermo- 
stat. This conclusion qualitatively agrees with some recent simulations of dense gases [Puglisi et 
\f~} ' al., 2007 J. Stat. Mech. P08016], although the deviations observed in computer simulations are 

t more important than those obtained here from the Enskog kinetic theory. Possible reasons for the 

quantitative discrepancies between theory and simulations are discussed. 
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I. INTRODUCTION 



The generalization of the fluctuation-response relation to non-equilibrium systems has received a considerable 
, attention in the past few years. In this context, granular matter can be considered as a good example of a system 
that inherently is in a non-equilibrium state. Granular systems are constituted by macroscopic grains that collide 
inelastically so that the total energy decreases with time. On the other hand, a non-equilibrium steady state (NESS) is 
reached when the system is heated by the action of an external driving force (thermostat) that does work to compensate 
' for the collisional loss of energy. In these conditions, some attempts to formulate a fluctuation-response theorem based 
on the introduction of an effective temperature have been carried out [H 0, 0, El • However, a complete analysis of 
the validity of the theorem requires the knowledge of the full dependence of the response and correlation functions on 
frequency iv Q . Given that this dependence is quite difficult to evaluate in general, the corresponding limit uj — > is 
usually considered. In this limit, the classical relation between the coefficients of diffusion D (autocorrelation function) 
. and mobility /x (linear response) is known as the Einstein relation. 
^ ' The Einstein relation for heated granular fluids has been widely analyzed recently. First, some computer simulation 
results for dilute systems have shown the validity of the Einstein relation (e = D/Tofj, = 1) in NESS when the 
1 temperature of the bath T is replaced by the temperature of the impurity Tq. This has an interesting consequence 
in the case of mixtures (where the different species have different temperatures @, B @) since a linear response 
. 1 experiment on a massive intruder or tracer particle to obtain a temperature measurement yields the temperature of 
the intruder and not the temperature of the surrounding gas. On the other hand, from an analytical point of view, 
kinetic theory calculations based on the Boltzmann equation have shown the violation of the Einstein relation (e ^ 1) 
in the free cooling case TO] as well as for driven granular gases [llj. These deviations are in general very small in 
the driven case (less than 1% when the system is driven by a stochastic thermostat) and are related to non-Gaussian 
properties of the distribution function of the impurities. This is the reason why such deviations cannot be detected 
in computer simulations of very dilute gases. 

However, a recent computer simulation study of Puglisi et al. [HJ at high densities has provided evidence that the 
origin of the violation of the Einstein formula is mainly due to spatial and velocity correlations between the particles 
that are about to collide rather than the deviation from the Maxwell-Boltzmann statistics. These correlations increase 
as excluded volume and energy dissipation occurring in collisions are increased. The simulation results obtained by 
Puglisi et al. [13] motivate the present paper and, as in the case of a dilute gas , kinetic theory tools will be used to 
analyze the effect of density on the possible violation of the Einstein relation. For a moderately dense gas, the Enskog 
kinetic equation for inelastic hard spheres [131 ] can be considered as an accurate and practical generalization of the 
Boltzmann equation. As in the case of elastic collisions, the Enskog equation takes into account spatial correlations 
through the pair correlation function but neglects velocity correlations (molecular chaos assumption) [Til ]. Although 



'Electronic address: vicenteg@unex.es URL: http://www.unex.es/eweb/fisteor/vicente/ 



2 



the latter assumption has been shown to fail for inelastic collisions as the density increases [TH, [ljl QjJ, there is 
substantial evidence in the literature for the validity of the Enskog theory for densities outside the Boltzmann limit 
(moderate densities) and values of dissipation beyond the quasielastic limit. This evidence is supported by the good 
agreement found at the level of macrosco p ic p roperties (such as transport coefficients) between the Enskog theory 
[H GJ, HI and simulation [HI [H, H, H S| and experimental [2g, l27| results. In this context, one can conclude 
that the Enskog equation provides a unique basis for the description of dynamics across a wide range of densities, 
length scales, and degrees of dissipation. No other theory with such generality exists. 



II. DESCRIPTION OF THE PROBLEM 



Let us consider a granular gas composed by smooth inelastic disks (d = 2) or spheres (d — 3) of mass m, diameter 
<7, and interparticle coefficient of restitution a in a homogeneous state. At moderate densities, we assume that the 
velocity distribution function /(v) obeys the Enskog kinetic equation [l3[. Due to dissipation in collisions, the gas 
cools down unless a mechanism of energy input is externally introduced to compensate for collisional cooling. In 
experiments the energy is typically injected through the boundaries yielding an inhomogeneous steady state. To avoid 
the complication of strong temperature heterogeneities, it is usual to consider the action of homogeneous external 
(driving) forces acting locally on each particle. These forces are called thermostats and depend on the state of the 
system. In this situation, the steady-state Enskog equation reads 

J7(v)=*7[v|/,/], (1) 

where J[v|/, /] is the inelastic Boltzmann collision operator, \ denotes the equilibrium configurational pair correlation 
function evaluated at contact, and T is an operator representing the effect of the external force. Two types of external 
forces (thermostats) are usually considered: (a) a deterministic force proportional to the particle velocity (Gaussian 
thermostat), and (b) a white noise external force (stochastic thermostat). The use of these kinds of thermostats 
has attracted the attention of granular community in the past years to study different problems. In the case of the 
Gaussian thermostat, T has the form [28L l29l. l30l| 

T!(y) = \ Q lk ' [V/(V)] ' (2) 

where £ is the cooling rate due to collisions. In the case of the stochastic thermostat, the operator T has the 
Fokker-Planck form [311 ] 

IT/ 9 x 2 



2 m \ov J 

The exact solution to the Enskog equation |T]) is not known, although a good approximation for / in the region of 
thermal velocities can be obtained from an expansion in Sonine polynomials. For practical purposes, one selects a 
finite number of terms in the expansion. In the leading order /(v) is given by 



/(v) - n-K-^v^e 



1 + - ( v* 4 - (d + 2)v 



(4) 



where v* = v/u t h, Vth — y/2T/m being the thermal velocity. Moreover, c is the fourth cumulant of the velocity 
distribution function / defined as 

r«*,_*ipi). (5) 



' did + 2) \4nT 2 

In the approximation ((4]) , cumulants of higher order have been neglected. Inserting Eq. ((4]) into the Enskog equation 
and neglecting nonlinear terms in c, one gets the following expression for the cooling rate £ [3lJ 



c = dr(f) xna Vth{1 ~ a ) V + 32°) ■ (6) 

The value of c depends on the thermostat used. In the case of the Gaussian thermostat, c is approximately given by 

m 

32(l-q)(l-2 Q 2 ) 
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FIG. 1: Plot of the temperature ratio Tq/T versus the coefficient of restitution a = ao for d — 3 in the case mo/m = oo/a = 0.5 
for the Gaussian thermostat (a) and the stochastic thermostat (b). The solid lines correspond to a dilute gas (0 = 0) while the 
dashed lines refer to a moderately dense gas (<f) — 0.2). 



while 

c(a) = 32(l-q)(l-2 Q 2 ) 
1 ' 73 + 56d-3(35 + 8d)a + 30a 2 (l-a) 1 1 

for the stochastic thermostat. It is interesting to remark that in the homogeneous problem the results obtained 
with the Gaussian thermostat are completely equivalent to those derived in the free cooling case when one scales the 
particle velocity with respect to the thermal velocity Vth In addition, although the expressions I©-© have been 
derived by neglecting nonlinear terms in the coefficient c, the estimates (J7J and ([8j present quite a good agreement 
with Monte Carlo simulations of the Boltzmann equation [30j, [32[ for moderate values of dissipation (say for instance, 
a ^ 0-5)- However, more recent results [33|, [HJ for the homogeneous (undriven) cooling state have shown that for 
very large inelasticity (a < 0.5), the higher-order cumulants may not be neglected since they can be of the same order 
of magnitude as c. The breakdown of the Sonine polynomial expansion is caused by the increasing impact of the 
overpopulated high-energy tail of the velocity distribution / [35| . 

We assume now that a few impurities or tracer particles of mass mo and diameter ao are added to the system. 
Given that their concentration is very small, the state of the gas is not affected by the presence of impurities. As 
a consequence, the velocity distribution function / of the gas still verifies the (homogeneous) Enskog equation (J]). 
Moreover, one can also neglect collisions among impurities themselves versus the impurity-gas collisions, which are 
characterized by the coefficient of restitution ao- Diffusion of impurities is generated by a weak concentration gradient 
Vno and/or a weak external field E (e.g. gravity or an electric field) acting only on the impurity particles. Under 
these conditions, the velocity distribution function /o(r, v, t) of impurities verifies the Enskog-Lorentz equation 



dtfo + v • V/ + — ■ ^-/ + Th = Xo J[v|/ 0l /], 
mo ov 



(9) 



where J [/o, /] is the (inelastic) Boltzmann-Lorentz collision operator and xo represents the equilibrium pair correlation 
function for impurity-fluid pairs at contact. Given that the gas is in a homogeneous state, it follows that xo is uniform. 
At a kinetic level, an interesting quantity is the partial temperature of impurities To. It is defined as 

^n To = J dv^V/ (v), (10) 

where uq is the number density of impurities. The corresponding cooling rate £o associated with the partial temper- 
ature Tq of impurities is defined as 

= ~T%T I dvm v 2 J[v\fo, /]. (11) 
dnolo J 



In the absence of diffusion (homogeneous steady state), Eq. © becomes 

^/o=XoJ[v|/ ,/]. 



(12) 
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FIG. 2: Plot of the coefficient Co versus the coefficient of restitution a = Qo for d — 3 in the case mo/m = uo/u = 0.5 for the 
Gaussian thermostat (a) and the stochastic thermostat (b). The solid lines correspond to a dilute gas ((j> = 0) while the dashed 
lines refer to a moderately dense gas (<j> = 0.2). 



This equation has been widely analyzed by using both types of thermostats [ll], [23| f° r bard spheres (d = 3). 
The results show that the temperatures of the gas (T) and impurities (T ) are clearly different and so the energy 
equipartition is broken down. In general, the temperature ratio 7 = Tq/T presents a complex dependence on the 
parameters of the problem. The condition for determining the temperature ratio 7 is different for each type of 
thermostat. In the case of the Gaussian thermostat, the temperature ratio is obtained by equating the cooling rates 

0.EJ 

C = Co, (13) 
while for the stochastic thermostat 7 is obtained from the condition [23| 



to mo 



(14) 



Requirements (|f 3[) and (|14|) lead to a different dependence of the temperature ratio Tq/T on the control parameters, 
namely, the mass ratio mo/m, the size ratio cto/ct, the coefficients of restitution a and ao, and the solid volume 
fraction 

na d . (15) 



2^dT(l) 

Apart from the temperature ratio, an interesting quantity is the fourth cumulant cq. It is defined as 



( ml f 4 d(d + 2) 



^0 = 37^(7^/ dvv%-^i). (16) 



As in the case of the coefficient c, the cumulant Co measures the deviation of fo from its Maxwellian form 

/o,M(v)^ 77 feV /2 ex P r^V (17) 



2T J "V 2T 

In order to determine the coefficients Co and cq one has to expand the velocity distribution function /o in terms of 
the orthogonal Sonine polynomials. As in the case of the granular gas distribution /, a good estimate of Co Emd cq 
can be obtained from the leading Sonine approach to fo: 



/o(v) -> 77 7T 1 v th 1 e 



(18) 



where 6 — m a T ' / 'mT is the mean square velocity of the gas particles relative to that of impurities. Expressions for (0 
and Co have been derived in Appendix [A] for an arbitrary number of dimensions d by considering only linear terms in c 
and cq. These expressions extend previous results derived in Ref. [Il| for hard spheres. Once Co & n d Cq are known, the 
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FIG. 3: (Color online) Plot of the Einstein ratio e — D/TofJ. versus the size ratio (Jq/o in the case of the Gaussian thermostat 
for d = 3 and a — ao = 0.5 when the impurities have the same mass density as the gas particles. Three different values of the 
solid volume fraction are considered: <f> = (solid line), <j> = 0.2 (dashed line), and <fi = 0.4 (dotted line). 



temperature ratio can be obtained from the constraints (fl~3"|) and (jT5J) for the Gaussian and stochastic thermostats, 
respectively. To get this explicit dependence, the form of the pair correlation functions x an d Xo m terms of the size 
ratio (Jo /a and the solid volume fraction (f> must be given. For a three-dimensional gas (d = 3), a good approximation 
for these functions is [33 
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where a = (a + oo)/2. For a two-dimensional gas (d — 2), \ and Xo are approximately given by [36 

X : 



(1 
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Xo 



9 a 
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(19) 



(20) 



(21) 



(22) 



16 a (1 - 0)2 

Obviously, x — Xo if <r = (To- Thus the temperature ratio and the kurtosis cq become independent of density for 
equal-size particles. The dependence of Tq/T on the (common) coefficient of restitution a = ao is illustrated in Fig. 
[T]for d = 3 in the case mo/m = <jq/(J = 0.5 and for two values of the solid volume fraction <f>. We consider the two 
types of thermostats discussed before. There is an evident breakdown of the energy equipartition in both thermostats, 
especially in the case of the stochastic driving force ([3]). However, the influence of density is more significant for the 
Gaussian thermostat than for the stochastic one. The dependence of Co on a is plotted in Fig. [2] for the same cases 
as considered in Fig. [TJ It is apparent that for both thermostats the value of Co is quite small for not too inelastic 
systems. This means that in this range of values of ao the distribution /o of the homogeneous state is quite close to a 
Maxwellian at the temperature of the impurity particle Tq. However, the magnitude of Co increases significantly as the 
dissipation increases, especially in the case of the Gaussian thermostat. This is a signal of the strong non-Maxwellian 
behavior of the reference homogeneous state of impurities for quite extreme values of inelasticity. As a consequence, 
cumulants of higher order than the fourth cumulant cq should be considered to assess the deviation of fo from its 
Maxwellian form /o,m for very small values of a and ao- Finally, with respect to the influence of density, we observe 
that it is more relevant for the Gaussian thermostat than for the stochastic thermostat, being practically negligible 
in the latter case. 



III. THE EINSTEIN RELATION 



The Einstein ratio e is defined as 



e = 



D 



(23) 
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FIG. 4: (Color online) Plot of the Einstein ratio e = D/Tofi versus the coefficient of restitution ot — ao in the case of the 
Gaussian thermostat for d = 2, mo/m = <tq/(T = 0.5 and three different values of the solid volume fraction: = (solid line), 
4> — 0.2 (dashed line), and <f> = 0.4 (dotted line). 



where D and /i are the diffusion and mobility coefficients, respectively. If the Einstein relation would hold, one would 
have e = 1. I want here to analyze the influence of density on e. The transport coefficients D and \x can be determined 
by solving the (inelastic) Enskog-Lorentz equation (J9j) by means of the Chapman- Enskog method [2j|. In the first 
order of the expansion, the current of impurities jo has the form [ll[ 

jo = -£»Vlnn + ^E. (24) 

Given that \ and Yp are uniform in this problem, it is evident that, when properly scaled, the previous solution 
obtained in Ref. [ll| for a dilute gas can be directly translated to the Enskog equation by making the changes 
J[fif] ~~ > xJ[fif] an d J[fo<£\—> XoJ[fo,f}- Technical details on the calculation of D and u by means of the 
Chapman-Enskog expansion [301 ] up to the second Sonine approximation can be found in Ref. for inelastic hard 
spheres (d — 3). The extension to an arbitrary number of dimensions is straighforward. Taking into account these 
results, the dependence of the Einstein ratio e on the parameter space of the problem can be obtained. In the case of 
the Gaussian thermostat, the result is 

e = 1 -y- — 37' (25) 

where the collision frequencies ^2 and 1/4 are explicitly given in Appendix [S] for an arbitrary number of dimensions d. 
The result for the case of the stochastic thermostat is 

e = l-™ (26) 
2 



It is clear that e becomes independent of the density when a — ao- As in the case of dilute gases [HI], Eqs. (j25|) 
and (|26p show that the violation of the Einstein relation in a heated moderately dense granular gas is basically due 
to the departures of fo from its Maxwellian form /o,m- At this level of approximation (second Sonine approximation 
to the coefficients D and fi), the deviation of fo from fo,M is only accounted for by the fourth cumulant cq. However, 
the coefficient cq depends on the solid volume fraction cj) through its dependence on the temperature ratio Tq/T. 
To assess the influence of density on the Einstein ratio e, Fig. [3] shows a plot of e versus the size ratio do/a in the 
case of the Gaussian thermostat for d — 3 and a = = 0.5 when the impurities have the same mass density 
as the gas particles [namely, mo/m = (<To/c) 3 ]- Three different values of the solid volume fraction <fi have been 
considered: = (dilute gas), = 0.2 (moderate dense gas), and = 0.4 (quite dense gas). The corresponding plot 
for the stochastic thermostat has not been included since the deviation of e from 1 is less than 1% for all the cases 
analyzed. It is apparent that the degree of violation of the Einstein relation is more important when the impurities 
are lighter and/or smaller than the gas particles, especially for high densities. To confirm this trend, the Einstein 
ratio has been plotted in Figs. H and [5] as a function of the (common) coefficient of restitution a = ao for d = 2, 
mo/m = Go/o = 0.5, and for the same values of as considered before. Figure [5] shows the results obtained by using 
the Gaussian thermostat and Fig. prefers to the results obtained for the stochastic thermostat. While e is close to 1 
in the case of the stochastic thermostat for all the densities considered, significant deviations form unity are observed 
for the Gaussian thermostat. In this latter case, it is apparent that the degree of violation of the Einstein formula 
increases with the volume fraction and the inelasticity. 



7 




0.98 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 5: (Color online) Plot of the Einstein ratio e — D/Tofi versus the coefficient of restitution a — ao in the case of the 
stochastic thermostat for d — 2, mo/m = ao/a = 0.5 and three different values of the solid volume fraction: <j> = (solid line), 
4> — 0.2 (dashed line), and (j> — 0.4 (dotted line). 
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FIG. 6: Plot of the ratio e(cj>)/e(0) versus the solid volume fraction <j) for a = ao = 0.5, mo/m — 2 and ao/a = 0.25 in the case 
of the Gaussian thermostat for spheres (d — 3) and disks (d = 2). 



This latter conclusion qualitatively agrees with the results obtained by Puglisi et at [12| from computer simulations 
since they observe a significant violation of the Einstein formula when excluded volume effects and dissipation are 
increased. However, at a quantitative level, the deviations observed by Puglisi et al. [13] are larger than those found 
here [see Fig. 3 of Ref. The quantitative disagreement between simulations [l2j and the Enskog theory results 

(l25f and (|26|) could be due to two different and independent reasons. First, as noted before in Section|TTl the expression 
of the fourth cumulant Co used here has been obtained by neglecting cumulants of higher order and considering only 
linear terms in c and co. On the other hand, according to the previous results [33'.. .34] obtained in the free cooling case 
for a monocomponent gas, the value of Co could change significantly for very strong dissipation when cumulants of 
higher order were taken into account in the form of the homogeneous distribution. In this context, the estimate of cq 
obtained in this paper could be not reliable for this range of values of inelasticity and so, the quantitative deviations 
of e from 1 observed in Figs. [3] and Q] for the Gaussian thermostat when a < 0.5 could be questionable. As a second 
reason of discrepancy between simulations and theory, one could argue that the velocity correlations (absent in the 
Enskog equation but present in computer simulations) play a more important role than spatial correlations (excluded 
volume effects) in the violation of the Einstein formula. In this case, one should correct the Enskog equation by 
incorporating recollision events ("ring" collisions) that take into account multiparticle collisions. 

Finally, since the results derived in this paper holds for a d-dimensional system, it is interesting to investigate 
the influence of dimensionality on the violation of the Einstein relation. To illustrate this effect, Fig. [6] shows the 
dependence of the ratio e(0)/e(O) on the solid volume fraction (f> for a = a — 0.5, mo/m = 2 and oq/o = 0.25 in the 
case of the Gaussian thermostat. I have considered the physical cases of hard spheres (d = 3) and hard disks [d = 2). 
Here, e(0) corresponds to the value of the Einstein ratio for a dilute gas. Although the qualitative dependence of the 
ratio e(0)/e(O) on <f> is quite similar in both systems, we observe that the violation of the Einstein ratio is stronger for 
d = 3 than for d = 2 for moderate densities. However, this trend changes as density becomes larger. 
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IV. CONCLUSIONS 



In this paper I have analyzed the validity of the Einstein relation e = D/Tq^j, = 1 for driven moderately d- 
dimensional dense granul ar g ases in the framework of the Enskog equation. This work extends a previous study 
carried out by the author [ll[ in the case of a dilute gas (</> = 0) of inelastic hard spheres (d = 3). To achieve a NESS, 
two types of thermostats (external forces) have been considered: (i) an "anti-drag" force proportional to the particle 
velocity (Gaussian force), and (b) a stochastic force, which give frequent kicks to each particle between collisions. 
The present work has been motivated by recent computer simulation results by Puglisi et al. [12j | where the spatial 
and velocity correlations between the particles have shown to be the most important ingredient in a strong violation 
of the Einstein relation. It is shown that e ^ 1, especially in the case of the Gaussian thermostat when the impurity 
is lighter and/or smaller than particles of the gas. As in the case of a dilute gas [TTJ] , the violation of the Einstein 
relation is connected to the strong non-Maxwellian behavior of the homogeneous velocity distribution function of 
impurities, which is mainly measured through its fourth cumulant c$. The results also show that the deviation of 
the Einstein ratio from 1 is more important as both the density and dissipation increase, which is consistent with 
the observations made by Puglisi et al. However, at a quantitative level, the deviations of the Einstein formula 

obtained here from the Enskog equation are smaller than those found in computer simulations. As discussed before, 
this quantitative disagreement between theory and simulation could be due to (i) the possible lack of convergence 
of the Sonine polynomial expansion for the reference Gaussian driven state, and/or (ii) the influence of velocity 
correlations which are absent in the Enskog description (molecular chaos hypothesis). With respect to the first source 
of discrepancy, one perhaps should include cumulants of higher order as well as nonlinear terms in c and Co to get an 
accurate estimate of the fourth cumulant cq for very small values of the coefficients of restitution a and «o- However, 
given that perhaps the absolute value of the higher order cumulants increases with inelasticity, the Sonine expansion 
could be not relevant in the sense that one would need an infinite number of Sonine coefficients to characterize the 
reference state. In this case, a possible alternative would be the use of the Direct Simulation Monte Carlo (DSMC) 
method (4p| to numerically solve the Enskog equation in the homogenous driven problem. Regarding the influence of 
velocity correlations, the inclusion of this new ingredient in the Enskog collision operator makes analytic calculations 
intractable since higher-order correlations must be included in the evaluation of the collision integrals. This contrasts 
with the explicit results reported in this paper, where the transport coefficients D and fi have been explicitly obtained 
in terms of the parameters of the system (masses, sizes and coefficients of restitution). 

Finally, it must be noted that the theoretical results derived here have been obtained by considering the second 
Sonine approximation to the Chapman-Enskog solution. Exact results can be obtained if one considers the inelastic 
Maxwell model (IMM) for a dilute gas. This model has been widely used by several authors as a toy model to 
characterize the influence of the inelasticity of collisions of the physical properties of the granular fluids. The fact 
the the collision rate for the IMM is velocity independent allows one to exactly compute the transport coefficients of 
the system. In particular, the coefficients D and /i have been evaluated [4l[ from the Chapman-Enskog method for 
undriven systems. The extension of such calculations to driven systems is straightforward. Thus, in the case of the 
Gaussian thermostat, one gets 

mo V 2 / m o V 2 / 



while 



D=*n, (28) 



in the case of the stochastic thermostat. Here, 



v D = -±— (l + a ), (=—(l-a 2 ), (29) 

d m + mo Id 

where w and Wq are effective collision frequencies of the model. According to Eqs. (|2~T|) and (f2"5|) , e = D/T /i = 1 
for both thermostats so that, the Einstein relation holds for the inelastic Maxwell model in any dimension. This 
conclusion agrees with previous independent results obtained for d = 1 0, |42j , d = 2 and d = 3 [|| . 
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APPENDIX A: EXPRESSIONS OF Co, Co, v 2 AND v A 

The explicit expressions of the partial cooling rate Coj the kurtosis cq and the collision frequencies v 2 and are 
displayed in this Appendix for an arbitrary number of dimensions d. In order to get these expressions, we consider 
the leading Sonine approximations ([U for the granular gas distribution / and (|18p for the impurity distribution f . 
The cooling rate Co can be obtained by following the same mathematical steps as those used before in previous papers 
[ill |3(| . The final expression can be written as 

Co = A 00 + Aoico + A 02 c, (Al) 



where 



47r (d-l)/2 /1 + 0\V2 

A °° = dirT*) XaUtT VthM \ o~ ) + 



M 

l-_(l + a o )(l + 0) 



(A2) 



Aoi - gdT jj; xon* d - 1 v th M^-^L (i + ao ) [2(3 + 49) - 3M(1 + a )(l + 0)] , (A3) 



.__- XoT ^-VM (^-J-J (l + ao)[2 + 3M(l + ao)(l + 0)]. (A4) 



Here, = m$T '/mTo and M = m/(m + mo). 

In order to get the coefficient Co, one substitutes Eqs. (|3]) and (TTSJ) into the Enskog-Lorentz equation (fT2|) . multiplies 
it by w 4 and integrates over the velocity. After some algebra and neglecting nonlinear terms in c and Co, the result in 
the case of the Gaussian thermostat is 

AoO + ^02 c + d(d+2) 

2 A 00 + A i + d{d+2) M Q V U 01 

while 

A o + a 02 c + 1 r£f 2 - ) M 2 9 2 {n 0Q + n 02 c) 

Co = — y (A6) 

A 01 + 3(3+2) M 9 "01 



for the stochastic thermostat. In Eqs. (|A5|) and (|A6|) . Mo = mo/(m + mo) and the quantities 

_(d-i)/2 (1 + 6*1 



fioi 



(d-l)/2 (l + 8)~ 1/2 

n °° = T7fj Xo^~^thM 2 M v ^ (l + o ) 

{/ ^ 2 I c J _|_ c 

-2 [d + 3 + (d + 2)9} + M (1 + a ) (1 + 9) 11 + d + , 
V a + 3 

-8Af 2 (1 + a ) 2 (1 + 6) 2 + 2Af 3 (1 + a ) 3 (1 + 6>) 3 } , (A7) 

j^Xon^-WM ^ Q J /2 (1 + a ) 

x {-2 [45 + 15d + (114 + 39rf)6» + (88 + 32d)6» 2 + (16 + 8d)6> 3 ] 
+3M (1 + a ) (1 + 8) [55 + 5d + 9(10 + d)9 + 4(8 + d)9 2 ] 

-24M 2 (1 + a ) 2 (1 + Of (5 + 46») + 30M 3 (1 + a ) 3 (1 + 9) 3 } , (A8) 

n (d~i)/2 (i + 9r 5/2 

X {2 [d - 1 + (d + 2)9} + 3M (1 + a ) (1 + 0) [d - 1 + (d + 2)9} 

-2AM 2 (1 + a ) 2 (1 + 8) 2 + 30M 3 (1 + a ) 3 (1 + 6») 3 } , (A9) 
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have been introduced. Equations (|A2|) — (|A4|) and (|A7[) - (|A9[) are consistent with the results @, [H[ obtained for hard 
spheres (d = 3). Once the coefficient Cq is given in terms of 7 = M /M9, the parameters of the mixture and the solid 
volume fraction, the temperature ratio 7 can be explicitly obtained by numerically solving the condition (fl3|) for the 
Gaussian thermostat or the condition (TH|) for the stochastic thermostat. 
Finally, the collision frequencies v 2 and ^4 are given by 



r (d-l)/2 



;/ 2 



-ti-l 



« tllX oA/(l + «o)^(l + er 1/2 , 



(A10) 



where 



;/ 4 



.4 



ttttTICF Vth YnM(l + an) I 

d(d + 2)r(f) x v o; Vi + ' 



.3/2 



A - (d + 2) 



1 



2Af 2 



2al - 



d + 3 



-OtQ 



+ d + l)[d + 5 + (d + 2)9} 



-D 



-M (1 + 6>) {/3<T 2 [(d + 5) + (d + 2)6>] [(11 + d)a 12 - 5d - 7] 
-(T^O + d(15 - 7a ) + rf 2 (l - a ) - 28a ] - [d + 2) 2 (1 - a Q )} 
+3(d + 3)/3 2 6>~ 2 [d + 5+ (d + 2)6>] + 2/?6 , ~ 1 [24+ lld + d 2 + (d + 2) 2 t 



-(d + 2)e- 1 [d + 3+(d- 



(d + 2)(l + 8)9- 2 [d + 3 + (d + 2)8], 



(All) 



(A12) 



B 



(d + 2)(1 + 2/3) + M(l + 0) {(d + 2)(1 - ao) - [(11 + d)«o - 5d - 71/36" 1 } 
+3(d + S)/? 2 ^" 1 + 2M 2 ^2a 2 - ^-^a + d + lj 0- l {\ + 6) 2 
-{d + 2)9- 1 {\ + 0). 



(A13) 



Here, (3 = M - M9 = M (l - 7 _1 ). For d = 3, Eqs. (|A10j) and (|A11|) coincide with those previously reported 
for hard spheres. 
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